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Abstract 



We investigate an extended chiral soliton model which includes tt,p,uj and a mesons 
as explicit degrees of freedom. The Lagrangian incorporates chiral symmetry and broken 
scale invariance. A scalar-isoscalar meson a is associated with a quarkonium dilaton field 
with a mass m a ~ 550 MeV. We show that the scalar field with anomalous dimension 
slightly changes the static and electromagnetic properties of the nucleon. In contrast, it 
plays a significant role in nucleon-nucleon dynamics and gives an opportunity to describe 
well the two-nucleon interaction. 
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1 Introduction 



In ref.JI]] Furnstahl, Tang and Serot (FTS) proposed a new model for nuclear matter and 
finite nuclei that realizes QCD symmetries such as chiral symmetry, broken scale invariance 
and the phenomenology of vector meson dominance. An important feature of this approach is 
the inclusion of light scalar degrees of freedom, which are given an anomalous scale dimension. 
The vacuum dynamics of QCD is constrained by the trace anomaly and related low-energy 
theorems of QCD. The scalar-isoscalar sector of the theory is divided into a low mass part that 
is adequately described by a scalar meson (quarkonium) with anomalous dimension and a high 
mass part (gluonium), that can be "integrated out", leading to various couplings among the 
remaining fields. The application of the model to the properties of nuclear matter as well as 
finite nuclei gave a satisfactory description. Further developments of the model ]2], ^ showed 
that the light scalar related to the trace anomaly can play a significant role not only in the 
description of bound nucleons but also in the description of heavy-ion collisions.lt was also 
shown that the anomalous cannot be due to an effect of nuclear density on the trace anomaly 
of QCD [|. 

Here a natural question arises: What is the role of this light quarkonium in the description 
of the properties of a single nucleon, when it is taken into account in topological nonlinear 
chiral soliton models, which are similar to the FTS effective Lagrangian on the single nucleon 
level? In the present paper we introduce a dilaton field with an anomalous dimension into the 
7rpo;-model M and investigate some properties of single nucleons which emerge as solitons in 
the sector with baryon number one (B = 1). 

It is well known that a scalar-isoscalar meson,0 the sigma, plays an important role in the 
nucleon-nucleon (NN) interaction especially within one-boson-exchange (OBE) models 0. 
We remark that the missing medium range attraction was a long standing puzzle in Skyrme 
like models. Lately it has been shown that M explicit inclusion of a scalar meson into the 
Skyrme model produces in a natural way the desired attraction. But these studies have two 
shortcomings. Firstly, the scalar meson used in such models |7j has nothing to do with OBE 
phenomenology, as it has a large mass and is identified with a gluonium state. Secondly, the 
Lagrangians used in refs. M do not include explicit omega mesons at all, which may "spoil" 
the mechanism of the attraction due to its strong repulsion.^ Therefore, it would be quite 
interesting to investigate the central part of the NN interaction when the light scalar (a) and 
cj-mesons are both taken into account explicitly. This is what is done here. In particular, 
it is important that if one is to properly describe the intermediate range attraction in the 
central NN interaction, the successful description of the single nucleon properties within the 
7ipu model should not be destroyed. We also note that in an soliton approach with explicit 
regulated two-pion loop graphs one is able to get the proper intermediate range attraction. In 
that case, however, one does not stay within a simple OBE approach any more (as done here) 
and also needs to calculate the modifications of the isovector two-pion exchange to the p and so 
on. For comparison, we mention that recent developments of the original Bonn OBE potential 
performed at Julich also include multi-meson exchanges leading to a renormalization of various 
interactions, couplings and cut-off parameters ||. The model we investigate is related closely 

4 Here and in what follows, we call it sigma meson for simplicity. Although we assume it to be some kind of 
quarkonium state, its precise dynamical nature is of no direct relevance for the following arguments. 

5 Of course, the stabilization of the soliton in refs.0 via higher derivative terms also leads to repulsion in the 
central NN interaction, but this is less easily be interpreted as single meson exchange. 
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to the OBE approximation of the NN force. 



2 The npuja model 



Including a a-meson by means of the scale invariance and trace anomaly of QCD into the 
7rpa;-model § can be done in terms of the following chiral Lagrangian of the coupled irpua 
system, 

C2 p -2a/d f2 -2a/d f2 -2a/d 

C = -5— d^a - TrLpIS - J -*— Tr[l, + r, + igfp^ + igurf+ 



3 1 f 2 m 2 e~ 3(T/d 
+~gu;,B» - z (u,ruT + p^p^) + Tr(U - V 



(1) 



d 2 Slml 
16 



1 - e~ ia/d I — + 1 



4a 



where the pion fields are parametrized in terms of U = exp (if -7?//^) and £ = y/U, left /right- 
handed currents are given by L p = U + d p U, Z M = £ + <9 M £, r p = £<9 M £ + , and the pertinent vector 
meson (p, uj) field strength tensors are p pv = d p p u — d u p p + g[p p x p u ] and = d^ojy — d^u^. 
Furthermore, the topological baryon number current is given by B^ = £ Ma ^ 7 TrL Q ,LgL 7 /(247r 2 ). 

In Eq. ([[]), So is the vacuum expectation value of scalar field in free space matter, f n is the 
pion decay constant (f n = 93 MeV) and g = g pwn is determined through the KSFR relation 
g = mjxplf^. The model assumes the masses of p and u mesons to be equal, m p = m w = m. 
The mass of the a is related to the gluon condensate in the usual way |], [1(J = 2^J~C^/ (dSo), 
where d is the scale dimension of scalar field (d > 1). Being "mapped" onto the states of a 
nucleon, the Lagrangian Eq.(P will be similar to the FTS effective Lagrangian. 

Nucleons arise as soliton solutions from the Lagrangian Eq.(|l|) in the sector with baryon 
number B — 1. To construct them one goes through a two step procedure. First, one finds 
the classical soliton which has neither good spin nor good isospin. Then an adiabatic rotation 
of the soliton is performed and it is quantized collectively^] The classical soliton follows from 
Eq.([T|) by virtue of a spherical symmetrical ansatze for the meson fields: 

G(r) 

U(f) = exp (ifrQ(r)), p a { = e iak r k , uJr) = uj(r)5 p0 , a(r) = a(r) . (2) 

gr 

In what follows we call 0(r), G(r), ou(r), and a(r) the pion-, p-, uj-, and cr-meson profile 
functions, respectively. The pertinent boundary conditions to ensure baryon number one and 
finite energy are, 9(0) = vr, G(0) = -2, u/(0) = cr'(0) = 0, 0(oo) = G(oo) = lj(oo) = cr(oo) = 0. 
To project out baryonic states of good spin and isospin, we perform a time-independent SU(2) 

6 We again refer the reader to ref. || for details. 
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rotation 

U(r,t) = A(t)U(r)A+(t), £(r,t) = A(t)£(r)A+(t) 
a(f,t) = air), u(r,t) = ^-[Kr\ 



(3) 



f ■ p {f,t) = -A{t)f- {KUr)+rK-fUr))A + {t) 
9 

f-p i {f,t)=A{t)f-p i {f)A+{t) 

with 2K the angular frequency of the spinning mode of soliton, if ■ K = A + A. This leads to 
the time-dependent Lagrange function 

C{t) = JdrC = -M H (B, G, w, a) + A(6, G, w, a, 0, £ 1; 6)Tr(ii + ) . (4) 

Minimizing the classical mass M#(0, G, u, a) leads to the coupled differential equations for 
0,G, io and a subject to the aforementioned boundary conditions. In the spirit of the large 
^-expansion, one then extremizes the moment of inertia A(G, G, uj, a, 0, £i, £2) which gives 
the coupled differential equations for £1,^2 and in the presence of the background profiles 
Q,G, uj and a. The pertinent boundary conditions are 0(0) = 0(oo) = 0, £i(0) = £1(00) = 0> 
^2(0) = £2(00) = 0, 2£i(0) + £2(0) = 2. The masses of nucleon and the mass of A , Ma, 
are then given by M N = M H + 3/8 A and M A = M H + 3/15A. 

The electromagnetic form factors are obtained in the usual way || are: 



3<? 

,-.<>\ 2nM N m 2 7 jAqr) _ 2(t/h 2 , 

G s M {q 2 = / — ^ r e 2 °' d r 2 dr, 

3gA J qr 

00 





,2 



^ > 



where s = sin(B),c = cos(G) and S2 = sin(G/2). The normalization is G'f(O) = G^(0) = 1/2. 
Similarly, meson-nucleon vertex form factors may be calculated | IT| . Of course, such strong 
interaction form factors are model-dependent quantities. In the soliton approach, however, they 
arise naturally as Fourier transfroms of the meson distribution within the extended nucleon. 
Stated differently, the soliton acts as a source of an extended meson cloud, which leads to a 
meson-nucleon interaction region of finite extension. In momentum space, this extension can 
be interpreted as a corresponding form factor. We mention that such a picture underlies the 
inclusion of strong form factors say in OBE models. Coming back to our approach, these form 
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factors are most easily evaluated in the Breit-frame: 



8rrM N f, 
3q 



oc 




77 



(q 2 + m 2 ) / ji(qr)sin(Q)r 



2 dr 



o 




— 20'c — B"rc + Q' 2 rs H h rm 2 s r 



r 



2 dr, 



oo 

a 



3o{qr) -& 



— - + m 4i p 1 

r 3 3r 



+ 



m 2 £, 





(6) 



where the "electric" and "magnetic" vector meson-nucleon form factors are connected to the 
Dirac Fi(t) and Pauli form factors i*2(0 through the following relations: G l E (t) = F[{t) + 
tFl(t)/AM 2 N , & M {t) = Fl(t) + F*{t), (i = p,u). 

3 Results and discussions. 

3.1 Static and electromagnetic properties of the nucleon. 

Using the formulas given above we have calculated static and electromagnetic properties 
of nucleon. As can be seen from Eq. QT|) , the Lagrangian has no free parameters in the 7rpu; 
sector. So, in actual calculations the parameters m^, m, f n are fixed at their emperical values, 
m w = 138 MeV, m = m p = m^ = 770 MeV, / w = 93 MeV, g = m/x/2/; = 5.85. In the cr-meson 
sector there are in general three free parameters: m a , Sq and the scale dimension d. The latter 
has been well studied for nuclear matter calculations |l], @ . In particular, it was shown that for 
d>2 the much debated Brown-Rho (BR) scaling may be recovered. Therefore, assuming that 
there is no dependence of d on the density, we shall use the best value d = 2.6 found in refs. fl], |2| . 
The values for So - were found to be So = 90.6 -r 95.6 MeV JT| . So we put So = f n = 93 
MeV. The mass of the a or equivalently the gluon condensate C g = m 2 (T d 2 Sl/'i = mid 2 f 2 /A is 
uncertain. We thus consider two cases: m a = 550 MeV and m a = 720 MeV in accordance with 
recent 7T7t phase shift analyses and with OBE values. We stress again that the precise nature 
of such a scalar-isoscalar field is not relevant here, only that it should not be a pure gluonium 
state. A summary of static nucleon properties obtained in both cases, i.e. with m a = 550 MeV 
and m a = 720 MeV, is given in Table [l]. One immediately observes that the nucleon mass is 
again overestimated. This may not be regarded as a deficiency, since it is known that quantum 
fluctuations tend to decrease the mass substantially. 

To estimate the influence of the cr-meson we also show the results given by minimal version 
of irpu model. As can be seen from Table [I], the inclusion of a light sigma meson into the basic 
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7rpo; model just slightly changes the nucleoli mass and its electromagnetic properties. This may 
be explained by the fact that the role of intermediate scalar-isoscalar meson in gamma-nucleon 
interactions is negligible. In contrast, the presence of the sigma-meson leads to an enhancement 
of axial coupling constant as it was first observed in ref . |I2| . Although, the physical mechanism 



of this change is not clear, it may be understood as mainly due to modification of meson profiles. 
This is in marked contrast to the inclusion of pion loop effects, which tend to lower the axial 
coupling even further In the present model qa = 0.88 for the Tipui and g& = 0.95 for the 



irpcua model, respectively. One may expect that an appropriate inclusion of a meson into a 
more complete version of the basic irpuj model might give the desired value Qa = 1.26. 

3.2 Meson— nucleon form factors and NN interaction 

One of the usual ways to calculate the meson-nucleon interaction potential within topolog- 
ical soliton models is the so-called product ansatz [[/[]. Within this approximation, the two- 
Skyrmion potential as a function of the relative angles of orientation between the Skyrmions 
has a compact form, and the extraction of the NN potential by projection onto asymptotic 
two-nucleon states is straightforward. This procedure gives only three nonvanishing channels: 
the central, spin-spin and tensor potential. At large and intermediate distances, the latter two 
compare well with e.g. the phenomenological Paris potential JL4). The major inadequacy found 



in such type of calculations is the lack of an intermediate range attraction in central potential. 
Although many remedies have been proposed, this result may not be genuine for Skyrme like 
models. In fact, the product ansatz, which is not a solution to the equations of motion, can 
only be considered accurate at large distances, and the failure of these calculations to reproduce 
the central range attraction may simply be the failure of the product approximation to provide 
an adequate approximation to the exact solution.^ Although the lack of central attraction may 
be recovered by the inclusion of a scalar-isoscalar meson, the inherent ansatz dependence of 
the trial configuration remains as a major shortcoming of product approximation [|l^ . 



On the other hand there is another natural way which was first used by Holzwarth and 
Machleidt [T7|]. They proposed to calculate Vnn within OBE model taking coupling constants 



and meson-nucleon form factors from a microscopical model such as the Cloudy Bag model or 
the Skyrme model. It was shown that the Skyrme form factor is a soft pion form factor that is 
compatible with the irN and NN systems. We shall use this strategy to investigate the NN 
potential within present model. 



The meson nucleon form factors given by well known procedure, proposed first by Cohen [TT 
are given in Eq.(|B|). Although they were derived in a microscopical and consequent way, these 
form factors could not be directly used in standard OBE schemes. The reason is that the OBE 
schemes || in momentum space use form factors defined for fields propagating on a flat metric, 
whereas the definition of form factors in Eq.(^|) involve a nontrivial metric. Hence, before using 
the latters in OBE scheme one should modify the procedure by redifining meson fields. The 
modification for pion nucleon form factors in 7rpu; model is clearly outlined in refs. |L7| . Now, 



7 An early study giving credit to this line of reasoning can be found in ref.flK 
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applying this procedure in the lagrangian in Eq. (|l|) we get the following ttNN form factor: 



3q 

OO 

nM N f n f ji{qr) 



(q 2 + m 2 ) / ji(qr)\J Mr{r) sm{Q)r 2 dr 



3 J qr 
o 



2F(r) 

-2F'(r) - rF"(r) + — + rm 2 F(r) 



r dr, 



(7) 



where 



M T = [1 + 2tan 2 (9/2)]e- 2,J / d , 
F(r) = 



3 + cos 2 (9(r)) - 4cos(0(r))e 



-a/d 



The influence of metric factor to the pion-nucleon form factor is illustrated in Fig. 1. It 
is seen that, without the inclusion of My the form factor is softer than in OBE models (the 
dashed line in Fig. 1), while its inclusion via Eqs. (0), (§) gives a behavior closer to OBE 
models. In fact, a monopole approximation at small q 2 = t of the normalised form factor 
G n (t)/G n (0) » A 2 /(A 2 - 1) gives = 860 MeV and A^ = 1100 MeV for M T = 1 and M T ^ 1 
respectively, compared to its emperical fit : A° BE = 1300 MeV (dotted line in Fig.l). Note, 
however, that our results for A^ are in line with recent coupled-channel calculations of the 
Julich group [0. There, a monopole form factor with A^ ~ 800 MeV is obtained. We do not 
want to stress here any qualitative comparison but rather like to point out that our approach 
also leads to cut-off values well below the ones obtained in OBE approaches. 

Introducing a flat metric requires a canonical form for the kinetic part of the lagrangian, 
which determines the dynamics of the field fluctation. The kinetic term of the scalar meson in 



Eq. © C 



kin 



S 2 e 



can be easily rewriten in a usual way: 

/Jin = Q^d^oj2 



by the following redifinition of the basic sigma field: a(r) = S d[l - e~ a ^' d }. Now the new 
field a may be identified with the real sigma field. Clearly this redifinition does not change the 
nucleons static properties given in Table |l|. Note also that, using the above redinition in the 
last term of Eq. (|IJ), one may easily conclude that m a = m a . The appropiate sigma-nucleon 
form factor is given by 



G°{~q 



4vr / jo(qr) 



2a' 



m 2 a 



r 2 dr, 



(9) 



and may be used in OBE models. We have not introduced any metric factors in the form 
factors of the heavier mesons since these should play a lesser role than in the case of the pion. 

For small values of the sqaured four-momentum transfer t each form factor can be parametrized 
in monopole form: Gi(t) = gi(A 2 — mf)/(A 2 — m 2 ) {i = n, p, lu, a) . We present in Table |2| the 
range parameters (cut-offs) and the coupling constants of the resulting meson-nucleon dy- 
namics. One can see that there is no interference between cx-meson nucleon and e. g. the 
pion-nucleon coupling constants. In other words, the inclusion of cr-meson does not signif- 
icantly affect meson-nucleon form factors that had been given by the irpu model. As it is 
seen from Table the values for meson-nucleon coupling constants are close to their emperical 
values (in some cases obtained by OBE model fits). This is one of the main advantages of the 
inclusion of a scalar-isoscalar meson as done in the present approach. 

The corollary of the present model is that it gives significant information on the u-nucleon 
interaction. As it is seen from the Table ||], the value for g a NN and the cut-off parameter of 
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sigma -nucleon vertex A a are smaller than their OBE prediction A° BE pa 1300-H2000 MeV. This 
contrast is evidently seen from the Fig. 2., where G a (t)/G a (0) for two cases : m a = 720MeV 
and m a = 550MeV is presented with the solid and dashed lines respectively. The band enclosed 
by the dotted lines refers to the OBE monople form factor with A° BE = 1300 . . . 2000 MeV. 
One can conclude that the present model gives a softer aNN form factor than it obtained 
by OBE. As it had been noticed before, the t-plane for each form factor has a cut along the 
positive real axis extending from t = to to oo. The cut for the cr-nucleon vertex function 
starts at to = 4m 2 reflecting the kinematical threshold for the a — > tttt channel. More precisely 
this result follows from the asymptotic behavior of the meson profiles: For r — > oo, we have 
O(r) ~ exp(— m 7r r)/m 7r r, a{r) ~ 2 (r), which are derived from the equations of motion. 

Once the vertex function of the corresponding meson-nucleon interaction has been found, its 
appropriate contribution to the NN interaction may be easily calculated by using well known 
techniques from OBE. The detailed formulas are given elsewhere || [T^|. In particular, the 
contribution of the cr-meson exchange to the central potential is given by 

CT 

The central NN potential in the T = 0, S = 1 state (the deuteron state) is presented in Fig. 3 
in comparison with Paris potential. Our prediction is in good agreement with the emperical 
one. Note that the desired attraction in the central Vnn has before been obtained in the irpuj 



model by means of two-meson exchange [I5|. 



In conclusion it should be noted that we do not intend (expect) to describe (cover) all NN 
phase shifts staying only in the framework of the present model. Besides other mesons, which 
are usually included in OBE picture, the full model should also take in account e.g. NAp 
couplings. In addition the 2n exchange and its strong mixing with a meson exchange (see e.g. 
ref. ||20|| ) should be considered. Another reason possible which limits the accuracy of iViV phase 



shifts in the present model is that the cxNN coupling is not sensitive the mass of sigma (see 
Table |2|) as it is in the OBE phenomenology. In fact, even when the 27r exchange is disregarded 
, the pure OBE model has to consider two types of sigmas with nearly the same masses but 
with quite different coupling constants. So, we refrain from performing direct calculations 
of NN phase shifts in the present model. Instead, we point out that the meson-nucleon form 
factors found in the present model could be useful in a wider context of calculations of nucleon- 
nucleon observables (phase shifts, deuteron properties etc) and may give more information on 
meson-nucleon and nucleon-nucleon dynamics. 

To summarize, we have developed a topological chiral soliton model with an explicit light 
scalar-isoscalar meson field, which plays a central role in nuclear physics, based on the chiral 
symmetry and broken scale invariance of QCD. We have shown that for the single nucleon 
properties, the successfull description of the electromagnetic observables of the irpu model is 
not modified and even the value for the axial-vector coupling is somewhat improved. In the 
two-nucleon sector, this extended npua Lagrangian leads to the correct intermediate range 
attraction in the central potential and a soft aNN formfactor for both values of sigma meson 
mass m a = 550 MeV and m a = 720 MeV. 

Acknowledgements: We would like to thank R. Machleidt for providing us with computer 
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Akademischer Austauschdienst (DAAD). 
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Table 1: Baryon properties in the npco and npua models 





TlpUl 


ixpoja 


TipUJG 


Exp. 


m a [MeV] 


— 


550 


720 




Cy 4 [MeV] 


— 


258 


295 


300^400 


5 1 / 4 [MeV] 


— 


119 


121 


— 


M N [MeV] 


1560 


1492 


1511 


939 


A [fm] 


0.88 


0.88 


0.88 


— 


M A - M N [MeV] 


344 


350 


350 


293 


th = {t 2 b ) 1/2 [fm] 


0.5 


0.5 


0.5 


-0.5 


(rlVj 2 [fm] 


0.92 


0.94 


0.94 


0.86± 0.01 


(r|)n [fm 2 ] 


-0.20 


-0.16 


-0.16 


-0.119±0.004 




0.84 


0.85 


0.85 


0.86± 0.06 


(r 2 M )l/ 2 M 


0.85 


0.85 


0.85 


0.88± 0.07 


//p [n.m.] 


3.34 


3.33 


3.33 


2.79 


fi n [n.m.] 


-2.58 


-2.53 


-2.53 


-1.91 




1.29 


1.30 


1.30 


1.46 


fi'A 


0.88 


0.95 


0.95 


1.26±0.006 


(ri)V^ [ fm] 


0.63 


0.66 


0.66 


0.65± 0.07 
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Table 2: Meson-nucleon coupling constants and cut-off parameters of meson-nucleon form 
factors. The A$(i = 7i,p,u)a) are cutoff parameters in equivalent monopole fits 1/(1 — £/A|) 
to the normalized form factors GWatW/GWaKO) around t — 0. The empirical values are from 
QBE potential fit 





TipijJ 


7rpu;o-(550) 


7rptU(T(720) 


OBE/Emp. 


GnNN(0) 


14.74 


13.97 


14.17 


13.53 


Go- AT AT (0) 


— 


6.2 


6.19 


9.1 (12.41) 


Ff(0) 


2.55 


2.76 


2.68 


2.24 




14.33 


15.01 


14.67 


13.7 




5.6 


5.43 


5.47 


6.1 


F?(0) 


8.78 


10.73 
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Fig. 1 . 




Figure 1: The normalized tcNN form factor in the npua model (m a = 720 MeV). 
The solid line represents the form factor when the metric factor is included (Eq. 
(0)), while the dashed line gives the result with no metric factor as in Eq. (^). The 
dotted line is a monopole form factor with = 1300 MeV. 
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Fig. 2. 
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Figure 2: The sigma-nucleon form factor G a NN(q 2 ) /GWjv(O)- The dashed and 
solid lines are for m a = 550 MeV and m a = 720 MeV, respectively. Typical OBE 
monopole fits with A = 1.3 ... 2 GeV are shown by the band enclosed by the dotted 
lines. 
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Fig. 3. 




Figure 3: The central potential in the S — 1, T — state for irpuj and irpua models 
(dashed and solid lines, respectively). No contribution from two-meson exchange 
has been taken into account. The dotted line corresponds to the Paris potential |n 
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